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We study the bifurcations and the chaotic behaviour of a periodically forced double- well Duffing 
oscillator coupled to a single-well Duffing oscillator. Using the amplitude and the frequency of the 
driving force as control parameters, we show that our model presents phenomena which were not 
observed in coupled periodically forced Duffing oscillators with identical potentials. In the regime of 
relatively weak coupling, bubbles of bifurcations and chains of symmetry-breaking are identified. For 
much stronger couplings, Hopf bifurcations born from orbits of higher periodicity and supercritical 
Neimark bifurcations emerge. Moreover, tori-breakdown route to a strange non-chaotic attractor is 
also another highlight of features found in this model. 
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I. INTRODUCTION 

The dynamics of coupled nonlinear oscillators has at- 
tracted considerable attention in recent years because 
they arise in many branches of science. Coupled oscilla- 
tors are used in the modeling of many physical, chemical, 
biological and physiological systems such as coupled p-n 
junctions charge-density waves 0, chemical-reaction 
systems [3|7 and biological-oscillation systems . They 
exhibit varieties of bifurcations and chaos [5-12], pattern 
formation |13| . synchronization [7,14-20] and so on. 

Bifurcation analysis is a useful and widely studied sub- 
field of dynamical systems. The observation of the bifur- 
cation scenario allows one to draw qualitative and quan- 
titative conclusions about the structure and dynamics of 
the systems theory. Numerous theoretical and numerical 
studies in this direction have been carried out for sev- 
eral specific problems [6-12]. Here, we have referred to 
only recent studies relevant with coupled nonlinear sys- 
tems. Of particular interest in this study are coupled 
Duffing oscillators - strictly dissipative nonlinear oscilla- 
tors which model the motion of various physical systems 
such as a pendulum, an electrical circuit or a Josephson- 
junction to mention but a few. Two or more identical 
or non-identical oscillators may be coupled in different 
ways. Thus, the type of the oscillators and couplings 
as well as the external perturbation have a significant 
influence on the dynamics of the coupled system. For 
instance, Kozlowski et al. @ have shown that the global 
pattern of bifurcation curves in the parameter space of 
two coupled identical periodically driven single- well Duff- 
ing oscillators consists of repeated subpatterns of Hopf 
bifurcations - a scenario that differs from that of a single 
Duffing oscillator. 



Very recently, Kenfack studied the model already 
considered in [6j but using identical double-well poten- 
tials. That model presented many striking departures 
from the behaviour of coupled single-well Duffing oscilla- 
tors. Scenarios such as multiple period-doubling of both 
types, symmetry-breaking, sudden chaos and a great 
abundance of Hopf bifurcations were found in that model, 
in addition to the well-known routes to chaos in a one- 
dimensional Duffing oscillator fTo| ]. 

The present study derives its motivation from refer- 
ences 0.0 . Here, we present a model consisting of a pe- 
riodically driven double- well Duffing oscillator coupled to 
a single-well Duffing oscillator. We numerically analyse 
varieties of bifurcations with special emphasis on those 
which are typical of the model. This is explored in dif- 
ferent coupling regimes as a function of the two control 
parameters of the system namely, the amplitude / and 
the angular frequency uj of the driving force. 

The paper is organized as follows: Section II describes 
the system under consideration. In section III, we present 
several numerical results of bifurcation structures and 
conclude the paper in section IV. 



II. THE MODEL 

The Duffing oscillator |2l| that we treat here is a 
well-known model of nonlinear oscillator. It is governed 
by the following dimensionless second-order differential 
equation 



d 2 x dx 
~d¥ + ~dt 



dV{x) 
dx 



/cos(wt), 



(1) 
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where x and b stand for the position coordinate of a par- 
ticle and the damping parameter, respectively. The right 
hand side of eq. (Q), represents the driving force at time t, 
with the amplitude / and the angular frequency lu. The 
oscillator belongs to the category of three-dimensional 
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dynamical systems (x, dx/dt, t). The system is a 
generalization of the classic Duffing oscillator equation 
and can be considered in three main physical situations, 
wherein the dimensionless potential 



2 4 

V{x)=a\+(3^, 



(2) 



is a (i) single-well (a > 0,(3 > 0), (ii) double-well 
(a < 0,(3 > 0) or (iii) double-hump (a > 0,0 < 0) 
potential. Each of the above three cases has become a 
classic central model describing inherently nonlinear phe- 
nomena exhibiting rich and baffling varieties of regular 
and chaotic motions. 

When two of such system eq.( |2J) interact with each 
other through a specific coupling, the dynamics is ex- 
pected to be even richer and more attractive. The cou- 
pling employed in the present paper is a linear feed-back 
which can be interpreted as a perturbation of each ocil- 
lator by a signal proportional to the differences of their 
positions. The potential governing such a coupled system 
reads 



V[x,y) = -x 2 



< [ +y2/ 2 + ?2/ 4 + |(--2/) 2 , (3) 



where k is the coupling parameter. Depending on the 
values taken by the parameters a, (3, a 1 and (3 , one dis- 
tinguishes three categories of globally bounded coupled 
Duffing oscillators; namely the single-single wells, the 
single-double wells and the double-double wells. Fig. 
displays a prototypical contour plot of the potential © 
for a weak coupling strength k = 0.1. Grey dots denote 
unstable fixed points (saddles), while dark ones repre- 
sent stable fixed points of the system. Clearly the single- 
double wells (Fig. nas an unstable fixed point and 
two stable fixed points; the single-single wells (Fig. ^b) 
shows only one stable fixed point; however the double- 
double wells (Fig.njc) possesses five unstable fixed points 
and four stable points. When the coupling strenght k in- 
creases, the positions of these fixed points are modified 
accordingly and provide a variety of interesting phenom- 
ena which we will study in the subsequent sections. From 
cc f(Gl? the equations of motion of the driven, coupled 
double-single wells (a = — a' > 0,(3 — (3' > 0) Duffing 
oscillators corresponding to (Fig. Ufa)) are derived and 
can be expressed as 



d x 
~d¥ 
fy 
dt 2 



dx 

-6— + ax — (3x 3 + k(y — x) + / cos(wi) 
dy 

-b— - ay - (3y 3 - k(y - x). 
dt 



(4) 



In the special case where k = 0, it is obvious that sys- 
tem Q reduces to two independent subsystems: the first 
being a periodically forced double-well Duffing oscillator 
with state variables (x,dx/dt), while the second is the 
unforced single-well Duffing oscillator having the state 
variables (y,dy/dt). The extended phase space of our 
model is five-dimensional (R 4 xS 1 ), wherein any element 
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FIG. 1: Contour plots of coupled Duffing oscillators for k = 
0.1 : (a) double-single wells, (b) single-single wells, (c) double- 
double wells. 



of the state space is denoted by (x, dx/dt, y, dy/dt, 9); S 1 
being the unit circle containing the phase angle, 9 = cut. 
Thus, in our numerical study, we visualize the attractors 
in the subspace along with their bifurcations by explor- 
ing the dynamics in the Poincare cross section defined by 



= {(x, dx/dt, y, dy/dt, 9 = 9 ) e M 4 x S 1 } 

where #o is a constant determining the location of 
the Poincare cross section on which the coordinates 
(x, dx/dt, y, dy/dt) = (x\,vi, X2,V2) of the attractors are 
expressed. It is worth noticing that time-reversal sym- 
metry is broken. Due to the symmetry of the potential 
V(x,y) = V(— x, —y), the system is invariant under the 
following transformation 

S : (xi,v 1 ,x 2 ,v 2 ,t) -> (-xi, -v-l, -x 2 , -v 2 ,t + T/2) 
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where T = 2tt/uj is the period of the driving force. There- 
fore in addition to saddle nodes (sn), period doubling 
(pd) and Hopf (H), one may expect symmetry breaking 
(sb). This can be evidenced with the eigenvalues spectra 
obtained from a simple linear stability analysis. Addi- 
tional features may arise not only because of the driving 
force but also because the two coupled oscillators are sub- 
ject to different potentials. 
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FIG. 2: (a) Bifurcation diagram for / = 4.0 and k = 0.1 
showing a sequence of reverse pd route to chaos. Periodic 
windows sandwitched by chaotic domains are also visible, (b) 
Lyapunov spectrum corresponding to (a). 



A. Weak coupling (0 < k < 1) 

At weak coupling strength, we find that the bifurca- 
tion structures are essentially similar to those commonly 
found in one-dimensional Duffing oscillators. For k = 0.1, 
we computed bifurcation diagrams in a large range of 
/. In Fig. |3 for instance, we set / = 4.0 and plot a 
typical bifurcation diagram in this regime, together with 
its Lyapunov spectrum X ma x- In the periodic window 
0.3375 < uj < 0.34, a sequence of reverse period-doubling 
(pd) bifurcations yields a period-5 attractor which is later 
destroyed at to — 0.34 in a crisis event, say sudden chaos. 
After the large chaotic domain intermingled with win- 
dows made up of periodic orbits, another reverse pd cas- 
cade occurs around 0.352 < uj < 0.355. Again, the cas- 
cade generates a period-4 attractor which is subsequently 
destroyed in another crisis event. This bifurcation sce- 
nario is in fact the common feature of the system as uj is 
further increased. However, for uj > 0.94, we find peri- 
odic orbits dominating the system behaviour. The tran- 
sitions are well characterized by the Lyapunov spectrum 
shown in Fig. |2Ib). 
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FIG. 3: Bifurcation diagram for k = 0.1 and / = 4.2. Bub- 
bles of bifurcation (at u) ~ 0.16) are sandwiched by period 
III. BIFURCATION STRUCTURES doubling (pd), born from symmetry-breaking (sb). 



In the numerical results that follow, we investigate 
the dependence of the system behaviour at given driving 
strength / and coupling strength fc, for varying angular 
frequency uj. The bifurcation diagrams show the projec- 
tion of the attractors in the Poincare section onto one of 
the system coordinates versus the control parameter. To 
gain further insight about the dynamics of the system 
under investigation, we compute the Poincare sections 
and Lyapunov spectra X m ax- These results are obtained 
solving eq.(EJ with the help of the standard fourth-order 
Runge Kutta algorithm. Starting from the initial condi- 
tion (xi, X2, V2) — (0.0, 0.5, 0.0, 0.5), the system is nu- 
merically integrated for 100 periods of the driving force 
until the transient has died out. To ascertain periodic, 
quasiperodic and chaotic trajectories, the system is fur- 
ther integrated for 180 periods. Next our results are pre- 
sented in weak, moderate and strong coupling regimes. 



As the driving-force amplitude / increases, different 
transition processes can be captured besides the pd se- 
quence already observed. The bifurcation diagram of 
Fig. plotted for / = 4.2, shows the occurrence of the 
attractor bubblings (bubbles) sandwiched by symmetry- 
breaking (sb). The sequence observed, consisting of a 
symmetry breaking, a period-doubling, bubbles, a re- 
versed period doubling and a symmetry breaking bifurca- 
tions, can simply be denoted by -sb-pd-(bubbles)-pd-sb-. 
This means at the first sb point, a symmetrical periodic 
orbit splits into two coexisting asymmetrical periodic or- 
bits (one being the mirror image of the other); each of 
the two asymmetrical periodic orbits give rise to another 
asymmetrical orbit with double period (pd) followed by 
bubbles. The bubbles in return generate asymmetrical 
periodic orbits which join by a pair (reversed pd) to 
form a symmetrical orbit at the second sb point. The 
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sequence -sb-pd-(bubbles)-pd-sb- which is also observed 
for / > 4.2, was not reported in previous literature for 
coupled Duffing oscillators. However, we found a similar 
sequence in a recent study for coupled ratchets exhibiting 
synchronized dynamics pj • Subsequently, we will refer to 
the integer n as the period of an orbit. 
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FIG. 4: Bifurcation diagram showing varieties of Hopf Bifur- 
cations for k — 1.0: (a) n = 2, / = 1 and (c) n = 6, / = 15.7; 
(b) Lyapunov spectrum characterizing (a). Here n represents 
the period of the periodic orbit. 



B. Moderate coupling (k — 1) 
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FIG. 5: Poincare sections (xi,Vi) illustrating the supercriti- 
cal Neimark bifurcation for k = 1 and / = 1: (a) u> = 0.38005, 
(b) uj = 0.38009 and (c) w = 0.3810. 



uj, the system moves outward from that unstable fixed 
point, passing through states like the one presented in 
Fig. Elb) for to = 0.38009, towards the attracting invari- 
ant closed curve (torus) shown in Fig.[5{c) for uj — 0.381. 
This transition is termed secondary Hopf bifuraction or 
Neimark Bifurcation [22ll23| . Moreover, it is a supercrit- 
ical Neimark bifurcation since the system moves outward 
from near an unstable fixed point towards the attracting 
invariant closed curve. In general, the bifurcated solution 
can be stable, say supercritical or subcritical otherwise. 



Setting k = 1, we find (beside pd and sb) an abun- 
dance of Hopf bifurcations. In the earlier work 6, 9], no 
periodic orbits with periods lager than n = 3 were found 
to undergo Hopf bifurcations. However in the present 
model, this number is much larger and can reach n = 6. 
In Fig. 01 we plot two bifurcation diagrams for (a) / = 1 
and (c) / = 15.7 in order to illustrate these observations. 
Hopf bifurcations with orbits of periods (a) n = 2 and (c) 
n — 6 can clearly be seen. Different transitions observed 
in (a) are well characterized by its Lyapunov spectrum 
displayed in (b). By means of Poincare cross-sections, 
we further zoom structures in this regime where Hopf 
bifurcations are predominant. Reported in Fig. [5] is an- 
other bifurcation scenario found in our model, namely the 
Neimark bifurcation. Fig.^a) displays an unstable fixed 
point spiraling out for uj = 0.38005. Further increasing 



C. Strong coupling (k = 5) 

Considering k = 5 as used in ref. 0, Q, and fol- 
low uj values, we find that periodic orbits dominate for 
large values of the driving amplitude (/ > 20). How- 
ever in Fig. phenomena such as period-doubling (pd) , 
symmetry-breaking (sb), saddle nodes (sn), resonance 
(R), bubbles and sudden chaos can be found in the fre- 
quency range 0.1 < uj < 0.4, for moderate values of 
/ < 20. Remarkably, the Lyapunov exponent charac- 
terizing these bifurcations exhibits large spikes at bi- 
furcation points. We have observed throughout a va- 
riety of Hopf bifurcations of higher period n < 6 (not 
shown), much like in the case k = 1. Interestingly, 
another transition found in this model is the passing 
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FIG. 6: Bifurcation diagrams for k = 5, (a) / = 20 and (b) 
/ = 4. Bubbles at around u = 0.26 and several regions of 
R, pd, sn, sb, and chaos are clearly seen; (c) the Lyapunov 
spectrum characterizing the bifurcations of (b). 

from low-dimensional to high-dimensional quasiperiodic 
orbits via tori-breakdown. The sequence of the trans- 
formation from tori (low-dimensional quasiperiodic or- 
bit) to a strange non-chaotic attractor (high-dimensional 
quasiperiodic orbit) are shown in Fig.[7| In Fig.[7{a), we 
plot a poincare section for u> = 1.38, showing a smooth 
quasiperiodic attractor visiting three folded tori in the 
(x 2 ,V2) sub-space. This folding state of tori is the on- 
set of tori-breakdown. Upon increasing the frequency to 
oj = 1.412, for instance, the smooth folded tori become 
completely broken, as can be seen in Fig. 0b). The cor- 
responding Lyapunov exponent, X max = —0.0174, shows 
that the attractor is not chaotic as the eye-test might 
tend to show; but rather it is a strange non-chaotic at- 
tractor. One would expect chaos from that transition as 
was reported by Kovlov [25| in a system of coupled Van 
der Pol oscillators. 



IV. CONCLUDING REMARKS 

In summary, we have introduced a model consisting 
of a periodically driven double-well Duffing oscillator lin- 
early coupled to a single-well Duffing oscillator. We have 
shown that the model admits a very complex dynamical 
structure that strongly depends on the coupling strength. 
In addition to the variety of phenomena earlier reported 
in the coupled Duffing oscillators with identical poten- 
tials 0, , our model exhibits several other interesting 
features which are basically attributed to the difference 
between the two oscillators's potentials. These are essen- 
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FIG. 7: Poincare sections (x-2,V2) showing the tori- 
breakdown route to a strange non-chaotic attractor; (a) at- 
tractor (folded tori) uo — 1.38, (b) a strange non-chaotic at- 
tractor (distroyed tori) u = 1.412, X„ la x = —0.0174. 



tially attractors bubbling (bubbles) , chains of symmetry- 
breaking bifurcations, Hopf bifurcations born form higher 
periodic orbits, supercritical Neimark bifurcations and 
the tori-breakdown route to a strange non-chaotic at- 
tractor. To visualise periodic, quasiperiodic and chaotic 
attractors of the system, Poincare cross-sections were 
utilized. Besides, the Lyapunov spectra used to char- 
acterise these transitions present spikes at the bifurca- 
tion points reflecting a sudden change in the system. In 
ratchet-like models, such transitions might lead to dra- 
matic changes in tr ansp ort properties of the system such 
as current reversals [24| . We hope these results would suf- 
ficiently complement previous ones and provide a more 
general overview, as far as bifurcations structures are 
concerned, in two coupled driven Duffing oscillators. Fi- 
nally, it would be interesting to examine broadly a situa- 
tion wherein the external forcing acts simultaneously on 
the two oscillators. 
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